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Baryon chiral perturbation theory (BChPT), as an effective field theory of low-energy quantum chromody- 
namics (QCD), has played and is still playing an important role in our understanding of non-perturbative strong 
interaction phenomena. In the past two decades, inspired by the rapid progress in lattice QCD simulations and 
the new experimental campaign to study the strangeness sector of low-energy QCD, many efforts have been 
made to develop a fully covariant BChPT and to test its validity in all scenarios. These new endeavours have not 
only deepened our understanding of some long-standing problems, such as the power-counting-breaking prob- 
lem and the convergence problem, but also resulted in theoretical tools that can be confidently applied to make 
robust predictions. Particularly, the manifestly covariant BChPT supplemented with the extended-on-mass-shell 
(EOMS) renormalization scheme has been shown to satisfy all analyticity and symmetry constraints and con- 
verge relatively faster compared to its non-relativistic and infrared counterparts. In this article, we provide a 
brief review of the fully covariant BChPT and its latest applications in the u, d, and s three-flavor sector. 

PACS numbers: 12.39.Fe Chiral Lagrangians, 12.38. Gc Lattice QCD calculations, 14.20.Gk Baryon resonances (S=C=B=0), 14.20. Jn 
Hyperons 
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I. INTRODUCTION 

Four fundamental interactions have so far been identified in nature: electromagnetism, weak interaction, strong interaction, 
and gravitation. Among them, the strong interaction is widely acknowledged as one of the most difficult to comprehend. In 
particular, the low-energy or non-perturbative strong-interaction phenomena have afflicted physicists for several decades. The 
non-abelian quantum field theory, generally accepted underlying the strong interaction, is the quantum chromodynamics (QCD) 
and was proposed about 30 years ago. However, two peculiar features of QCD makes its solution in the low-energy region 
difficult. First, the running coupling constant, a s , becomes large as the energy (or four momentum) transfer decreases-a phe- 
nomenon called the asymptotic freedom |H|2l. The large running coupling constant renders a perturbative treatment of QCD 
useless. Another closely related feature is color confinement. That is, although the degrees of freedom in QCD are quarks and 
gluons, which carry a quantum-number color, in nature only color-neutral objects, mesons and baryons, are observed. This leads 
to the complication that in QCD the underlying degrees of freedom and the observables are different. As a result, in studies of 
the low-energy strong interaction phenomena, despite the fact that there exists a fundamental theory, one has to turn to effective 
field theories (EFTs) or models, which are motivated by the various symmetries of QCD, or brute-force numerical methods, such 
as lattice quantum chromodynamics (LQCD) J^. 

One of such EFTs is the chiral perturbation theory (ChPT). In the massless limit, the QCD Lagrangian is invariant under the 
separate transformation of left and right handed quark fields, i.e., £q CD is invariant under SU(AT) l x SU(N)r (N = 2 or 3), 
which is termed as chiral symmetry. On the other hand, chiral symmetry is spontaneously broken, resulting in the appearance 
of eight massless Nambu-Goldstone bosons (NGBs). Experimentally the lowest-lying pseudoscalar octet can be identified as 
the corresponding NGBs. Furthermore, the masses of the u, d quarks and, to a lesser extent, that of the strange quark are small 
enough so that a perturbative treatment of their masses is appropriate 0|. 

It was Weinberg who first proposed that the soft-pion results of current algebra can be recast into the language of EFTs ijH . 
It was then systematically developed in the mesonic sector with two flavors |@] and three flavors 10]. Its application to the 
one-baryon sector came a few years later and met an unexpected difficulty, i.e., the combination of the power-counting (PC) rule 
of Weinberg and the modified minimal subtraction scheme (MS) yield the so-called power-counting-breaking (PCB) terms fl, 
which invalidates a systematic chiral expansion. To remove these PCB terms and restore the proper chiral PC rule, several 
approaches have been proposed. The first and most extensively studied is the heavy-baryon (HB) ChPT jgf]. It has been suc- 
cessfully applied to study a variety of physical phenomena, particularly in the u and d two-flavor sector (see, e.g., Ref. ifioll for 
a review of early applications). However, from the beginning, the drawbacks and limitations of the HBChPT are well known. 
Over the years, a number of relativistic formulations of BChPT were proposed by different authors llTTj - ll3ll . Among these, the 
infrared (IR) (H and extended-on-mass-shell (EOMS) lfl3ll schemes are the mostly studied. The infrared approach restores 
the PC rule in such a way that analyticity is broken at a scale of twice the baryon mass, which was originally believed to be of 
minor importance lfl2Tl . However, later studies found that the effects start to show up at much smaller meson masses, which are 
relevant for studies of the light-quark mass evolution of certain observables, such as the nucleon magnetic moment ifTill . and for 
studies performed in the three-flavor sector lll5ll . The EOMS approach, on the other hand, is free of such problems and appears 
to be a better formulation of covariant BChPT. Apart from the fact that a fully covariant formulation of BChPT is formally more 
appealing, it has been demonstrated that the EOMS approach converges relatively faster, a very important and practical feature 
for an EFT (see, e.g., Refs. QHH). 

In the past decade, lattice QCD has dev elop ed into an indispensable tool to study low-energy strong interaction phenomena 
(for some recent reviews, see, e.g., Refs. IU8U24II ). To reduce computing cost, most LQCD simulations employed larger than 
physical light-quark masses, finite volume and lattice spacing. As a result, to obtain the physical value of any observable 
simulated in the lattice, extrapolations to the physical world are necessary. ChPT provides an useful framework to perform such 
extrapolations and to estimate the induced uncertainties ll25ll . On the other hand, the LQCD quark-mass, volume, and lattice- 
spacing dependent results also offer a unique opportunity to help determine the many unknown low-energy constants (LECs) 
of ChPT 12611 . In many recent studies, it has been shown that the EOMS ChPT can provide a better description of the LQCD 
quark-mass dependent results than its non-relativistic counterparts lfl6[[l7ll . 

This paper is organized as follows. In Section II, we write down the relevant chiral Lagrangians needed for later discussions, 
introduce the PCB problem encountered in constructing a covariant BChPT, explain briefly the HB, IR, and EOMS formulations, 
and point out their advantages and limitations. In Section III, we show a few recent applications of the EOMS BChPT in the u, 
d and s three-flavor sector, including the magnetic moments of the octet and decuplet baryons, the masses and sigma terms of 
the octet baryons, the hyperon vector couplings, the recent developments in the u and d two-flavor sector, and the extension to 
the heavy-light system. Section IV contains a brief summary and outlook. 



II. COVARIANT BARYON CHIRAL PERTURBATION THEORY 

In this section, we briefly explain the HB, the IR and the EOMS formulations of baryon chiral perturbation theory (BChPT). 
This article is intended to be neither a pedagogical introduction to ChPT nor an exhaustive collection of all the recent works. For 
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both purposes, there exist excellent monographs and review articles flol, I27l433ll . Furthermore, we concentrate our discussions 
on the one-baryon sector in the three-flavor space. 

The chiral Lagrangians in the one-baryon sector can be generically written as 

C = C-MB + Cm, (l) 

where Cmb contains the part of the Lagrangians describing the interaction of the low-lying baryons with the NGBs while Cm 
describes the self-interaction of the NGBs. In general, the Cmb and Cm contain a series of terms ordered by the so-called chiral 
power counting (PC), i.e., 

Cmb = C M \ + C M \ + cfj B H , (2) 



Cm = C m ) +C m ) +C m ) + ---. (3) 

The chiral order, indicated by the superscript, is assigned in the following way: the mass of a NGB and a derivative on its field 
are counted as of 0(p). The mass of a baryon is counted as of O(p ) and the same is for the derivative on its field. However, 
f — tub, with P M the baryon four momentum, is counted as of 0(p). As a result, the baryon propagator j > J m is counted of 

as 0(p^ 1 ) and the meson propagator p , j 1 a as 0(p~ 2 ), with p denoting a generic small quantity, tub and itim the masses of 

the baryon and the meson. For a Feynman diagram consisting of L loops, Nb baryon propagators, Nm NGB propagators, and 
m k-th order vertices, its chiral order n is defined as 

n = AL - 2N B - 2N M + ^ mk. (4) 

With the chiral PC rule defined above and together with the MS scheme, this completes the definition of a PC rule for ChPT in 
the mesonic sector, i.e., in studies of the self-interaction of the NGBs. 

Once a matter field, such as the baryon field discussed here, is introduced, the above PC rule is violated. That is to say, 
lower-order analytical terms appear in a nominally higher- order calculation. For instance, in the calculation of the nucleon self- 
energy at 0(p 3 ), one finds (D(p ) and 0(p 2 ) terms. Although these terms do not change the underlying physics, they have two 
unwelcome consequences for an EFT. 1) They can completely change the natural values of the lower-order LECs. 2) They render 
an order-by-order analysis difficult. Both complicate a lot the study of baryon properties. In the following, we will explain in 
detail how one can recover the proper PC rule defined in Eq. (0} . 

A superficial difference between the Lagrangians (|2) and (|3) is that Lagrangian (01 contains only even-order terms while 
Lagrangian (fJJ has both even and odd order terms starting from order 1. It has a direct consequence: BChPT converges slower 
than meson ChPT. Furthermore, there are more LECs in the baryon ChPT than in meson ChPT. For instance, in the mesonic 
sector, there are two LECs at leading order (LO) and 12 LECs at next-to-leading order (NLO) 101:0 



4? = ^Tr[^C/(^C/)t] + ^Tr[ X U^ + U X % (5) 



C M ] = L x |Tr[^C/(^C/)t]} 2 + L 2 Tr[D f ,U(D„Uy]Tr[D»U(D»U^] 

+L 3 r Tr[D^U{D^UyD,U{D l/ Uy} + L 4 Tr[D AI C/(^C/)t]Tr[xC/ t + U X f ] 

+L 5 Tr[D^U(D»U)H X U 1 + U X ^)] + L 6 {Tr[ X tf + Ux^}} 2 

+L 7 {Tr[ X C/ t - UxHY + L 8 Tr[U X f U X 1 + X^W] 
-iL 9 Tr\f*D»U{iy'U? + fa(D»U?B"U\ + L w Tr[Uf^f^} 

+JTiTV[/*,/£' + f^JH + #2Tr[ XX t] } (6) 
where U = exp («-#-), D^U = d^U - ir^U + iUl^, f& = d^r v - d v r^ - i[r^,r„], /£, = cy„ - e?„Z M - i[Z M ,Z„], 



'F a J' A* — M M 1 ' j ^ 

X = 2i?o(s + ip), and $ collects the NGB fields 



/ -L ?? + 7T° V2tt+ V2K 



V2tt- -j= V - ir° V2K° I . (7) 



\ V2K- V2K° —fa 



1 The 0(p 6 ) Lagrangians can be found in Ref. |34|| . 



(a) 



FIG. 1. Feynman diagrams contributing to the nucleon mass at NLO BChPT. Diagram (a) represents tree level contributions, both at order 
and order 2, while diagram (b) represents contributions at NLO (order 3). The solid line denotes the nucleon and the dashed line denotes a 
pseudoscalar meson. The solid dot denotes an order 1 vertex coming from the Lagrangian of Eq. ((8). 

Presently, the LECs can be determined using either empirical inputs, estimates based on the resonance saturation assumption, or 
lattice QCD. See Ref. l35tl for a latest review on the present status in the mesonic sector. 

In the one-baryon sector, the LO Lagrangian has the form 

d$ B = B(i0 - m Q )B + -B^ 5 [u M B] + ^B 7 ^ 5 {u M B}, (8) 

which contains two LECs, D and F. At NLO, the number of LECs increases to 16 ll36Tl : 

C M \ = b D (B{ X+ ,B}) +b F (B[ X+ ,B}) +b (BB)(x+) 
+b 1 (B{u», [u^,B}}) + b 2 (B{u", {u M , B}}) 
+b 3 (B{u^[ Ut ,,B]}) +b i (BB)(u f *u ti ) 

+ib 5 [u" )7m D v B]]) - (bXK, K,7m£]]>) 

+ib 6 ((b[/,k, 7 ^b}]) - (bKk.K.t^I}}) 

+ib 7 ({B{u^{u v ,^D u B}}) - (bKk,!""^^}})) 

+ib s ((B^D„B) - (B*D vlll B)) {uV) +idi{B{[u»,u v ],a llv B}) 
+id 2 (B[[u»,u%a^B]) + idz(Bu> 1 )(u v a^B) 

+d 4 (B{f^,a^B})+d 5 (B[f^,a fw B}). (9) 

At the third order, the number of independent LECs increases to 78 l36l - l38ll . In general, the number of LECs increases rapidly 
with each increasing order-a feature of non-renormalizable EFTs. Fortunately, not all the LECs contribute to a particular 
physical process and therefore even at higher orders, BChPT retains its predictive power for certain observables, such as the 
hyperon vector form factors at zero four momentum transfer, where no LECs appear up to 0(p 4 ) 13911 . 



A. Power-counting-breaking and its restoration 

The chiral PC rule of Eq. (|4) is known to be violated in the one-baryon sector JD). In the following, we will use the nucleon 
mass as an example to demonstrate the PCB. Up to NNLO, the nucleon mass receives contributions from both tree and loop 
diagrams (see Fig. [TJ. Schematically, it reads as 



3(D + F) 2 

m N = m + 2(6 + 2b F )ml + \ ' loop^m^) 



m + cml + iff loop (3) (m ff ) (10) 
64^0 it 1 
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where loop^ 3 -* indicates the one-loop contribution from diagram (b) of Fig. Q] According to the chiral PC rule of Eq. (|4), diagram 
(b) counts as of order 3.0A direct calculation of the contribution of diagram (lb) to the nucleon mass yields 



loop 



(3) 



d 4 k 



jg(jg - f + m )# 



(2^) 4 (fc 2 - m\ + ie)((p - kY -ml + ie) 



(11) 



2ml 



logf 



2m ml 



log 



/'" 



1 

mo 



log 



2m; 



4m 2 



m; arccos 



\2m) 



The above result is divergent and therefore in the second line we have utilized the MS scheme to remove the divergence. Clearly, 
according to the chiral PC rule of Eq. (|4), the underlined terms are of order and 2, respectively, and therefore do not satisfy the 
specified chiral order of this diagram, C(p 3 ). These terms can be removed in a systematic way through a number of different 
approaches, such as the HB J^], the IR |[l2ll . and the EOMS lfl3tl schemes. The explicit expressions of the nucleon mass up to 
0(p 3 ) in these three approaches can be found in Ref. ll40ll . In the following, we briefly summarize the essential features of the 
HBChPT d, IR BChPT (H, and EOMS BChPT fill . 



B. The heavy baryon (HB) approach 

In the heavy -baryon (HB) approach H, one separates the nucleon four- momentum p into a large piece and a soft residual 
component 

p fi = mv fi + ^. (12) 

The velocity four-vector satisfies 

v 2 = 1, v° > 1, (13) 

which can be taken to be v* 1 = (1,0, 0, 0) for practical purposes. Note that with this choice of v 1 - 1 , the soft residual momentum 
satisfies = — fc 2 /(2m) = E—m<€i m. The essential idea is that although the zero-component of the nucleon four momentum 
is not small, the soft residual component k s is small compared to the nucleon mass. 

One can also decompose the relativistic nucleon field into two velocity-dependent fields 

1>(x)=e- iTm '-*(H v (x) + L v (x)), (14) 

where L v = e +tmv ' x P v+ ip and H v = e +unv ' x P v _ip are called light and heavy components, respectively, with the projection 
operators 

P v ± = ^(l + ^). (15) 

It can be easily shown that fL v = L v and fH v = —H v , From the relativistic Lagrangian of Eq. (|8), employing either the 
equation of motion technique or the path integral technique (see, e.g., Ref. 14111 ). one can obtain the corresponding LO HB 
Lagrangian 

= L v (iv ■ D + g A S v -u)L v , (16) 

where ~ indicates that the Lagrangians is in the HB formalism and Stf = i'ysC^Vp is the spin matrix. Note that the nucleon 
mass has disappeared from the LO Lagrangian and only appears in terms of higher orders as powers of 1/m. The HB propagator, 
derived from the above Lagrangian, has the following form 

G v (k)= f" + (17) 
v ■ k + i0 + 



The complete HB Lagrangians up to and including order p A can be found in Ref. (42 

Using the HB Lagrangians and propagators together with the dimensional regularization scheme, MS, yields the HBChPT. It 
has been widely applied to study various physical observables and in most cases turns out to be successful. For a comprehensive 
review of early applications, see, e.g., Ref. IU011 . 



2 We limit our discussion to the two flavor case and neglect the contribution of the virtual A(f 232). A complete SU(3) NNLO calculation in EOMS BChPT 
can be found in Ref. (3] . 
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On the other hand, it was noted in the very beginning that the 1 / m expansion can create analyticity problems under specific 
kinematics or, in other words, fails to converge in part of the low-energy region. For instance, the nucleon scalar form factor 
calculated in HBChPT explodes in the vicinity of t = 4m 2 . lllolfll. In addition, in HBChPT one has to keep track of various 
1/m corrections resulting from the non-relativistic reduction of the relativistic Lagrangians and, therefore, the number of terms 
in the effective Lagrangians increases more rapidly than in its covariant counterpart. 

To overcome the drawbacks of the HBChPT and derive a meaningful relativistic formulation of BChPT, many efforts have 
been taken in the past two decades. These efforts have deepened our understanding of the PCB problem and provided a number 
of useful solutions. Two most widely used formulations are the infrared (IR) lfl2ll and the extended-on-mass-shell (EOMS) lfl3ll 
schemes. 



C. The infrared (IR) scheme 

The IR scheme was introduced by Becher and Leutwyler lfl2ll . following the approach of Tang and Ellis ifTTll . It is based on 
the observation that in D dimensions, the infrared singular part of a loop diagram, which is free of PCB pieces, can be separated 
from the regular part, which is a polynomial of the momenta and quark masses. The regular part contains, in addition to the PCB 
terms, analytic higher-order pieces, which in principle can be absorbed by a redefinition of the LECs. To show the origin of the 
infrared singularities and how to remove the PCB terms in the IR scheme, we use the following integral as an example 

f d D q 1 

G = l J (2tt) d [(P -q) 2 -m 2 + ie] (q 2 - + ie) ' (18) 

where D denotes the number of space-time dimensions, mo and m w are the chiral limit nucleon mass and the lowest-order 
pion mass. This integration can be easily performed by use of the Feynman parameterization method. Introducing a Feynman 
parameter z and by use of the following identity, 



1 f 1 dz 



ab J [az + 6(1 — z)\" 



(19) 



one obtains 



f _d^q_ T dz 
l J (27tW ((P-zq^-M 2 ) 2 ' {M> 

where M 2 = P 2 z(z — 1) + m\z + m^(l — z) — ie. The integration over q can be easily performed, yielding 

G = r ( 2 - D / 2 ) jf dz(M 2 ) D/2 - 2 - (2D 

It was shown by Becher and Leutwyler that the integral G can be divided into two parts, the infrared singular part / and the 
remaining regular part R, G = I + R, defined as 

i r°° 
I = -J^yJj2 T ( 2 - D / 2 )J o dz(M 2 ) D/2 - 2 , (22) 

R = J^dJ2 r ( 2 - D l 2 ) ] dz ( M2 ) D ' 2 ~ 2 - C 23 ) 

It can be shown that for non-integer D the integral / is proportional to a non-integer power of the pion mass ~ M D ~ 3 and 
satisfies the PC rule of Eq. @. On the other hand, the remaining regular part R contains non-negative powers of the pion mass 
and violates the PC rule, for arbitrary D ITLTtl . The IR scheme then dictates that one takes into account only the infrared part and 
drops the regular part, with the argument that their contributions can be absorbed by the corresponding LECs. 
Calculating these three integrals explicitly in D dimensions and then taking the limit of D — > 4, one obtains 



„ 1 / 777, n \ 777?. , / 777^ \ 2m^ / 777?. / 777^ \ 

G = 77^ 1 - 2 + lo § [-I + "I lo § — + — 1 \ X ~ arccos [it- U ' (24) 



A finite-range regularization (FRR) scheme has been proposed to improve the convergence behaviour of HBChPT HHH. It has been rather successfully 



applied to study a variety of physical observables, particularly, in connection with LQCD simulations, their light-quark mass dependences (see, e.g., Refs 
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G 



IR 



16tt 2 



2m v 
m 



f to£ / m„ 

1 — ■ — 7T arccos 



Am\ 



2TOn 



m n my ^ / to*- 
2m 2m \ m 



(25) 



fi = TF~2 \ ~ 2 + lo S + — ff \/ 1_ T^ arccos 1 

I67H I V A* / m o TO o V 4toq \ 

where the conventional MS procedure has been applied to remove the divergent pieces. According to the chiral PC rule of 
Eq. (HJi, the two underlined terms in the above equations break the PC and therefore must be removed. In the IR scheme, not 
only these two terms, but also a series of higher-order terms are dropped. This can be clearly seen by expanding the regular part 

in TOjr, 

167T 2 ( ^ ' ( /i 2 ) TOq 6mg ) ^ ^ 

Furthermore, the infrared integral / is not analytical at to* = 2too because of the non-analyticity of the square root 

V^to^-to 2 . 

The infrared scheme has been reformulated by Schindler and collaborators Il54ll in a form analogous to the EOMS scheme and 
has been applied to study the nucleon mass up to 0(p 6 ) lf55l,l56tl . the axial, induced pseudoscalar, and pion-nucleon form factors 
up to 0(p 4 ) lf57tl . the electromagnetic form factors of the nucleon up to C(p 4 ) lf58ll . To overcome the non-analyticity induced by 
the original IR scheme and to establish direct connections between the values of the LECs determined in the HBChPT and those 
obtained in the co variant BChPT, Gail et al has introduced the so-called modified infrared renormalization scheme, IR f59ll . It 
has recently been applied to study the finite-volume corrections to the electromagnetic current of the nucleon ll60ll and the sigma 
terms and strangeness content of the octet baryons ll6lll . 
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D. The extended-on-mass shell (EOMS) scheme 



The essential idea of the EOMS approach is that one can calculate loop diagrams as usual in quantum field theory but then 
drop the PCB terms. Physically, it is equivalent to a redefinition of the existing LECs. This procedure is always valid since 
by construction ChPT contains all possible terms allowed by chiral symmetry, unitarity, and etc. For instance, in Eq. (TTOb the 
underlined PCB terms can be absorbed by too and c. It should be noted that this is different from the IR and HB procedures, 
because they remove not only the PCB terms but also a series of higher-order terms, which cannot be absorbed by the available 
LECs in question. 

Technically, one has a few different ways to implement the EOMS scheme. One can first calculate a Feynman diagram, obtain 
the analytical results and then subtract the PCB terms. A second alternative is to expand the integral in terms of some suitably 
defined small variables, perform the integration and obtain the PCB terms lfl3ll . One can also perform an expansion in terms of 
1/to and identify the relevant PCB terms. All three different implementations yield the same EOMS regularized result. In the 
case of the integral Eq. < fT~8T >, one obtains 

n 1 i m l i { m ir\ , 2m * L ™% ( 

Geoms = 2 \ — 2 lo S + \ 1 ~ arccos o — r ' ( 28) 

lo7r^ I TOg V Too / too y 4toq \ ztoo J I 

which has the same analyticity as the original MS result of Eq.([24l>. 



E. A brief summary of the different formulations of BChPT 

In Table U we summarize the main features of the four different formulations of BChPT. All of them, except the MS scheme, 
satisfy the PC rule of Eq. ©. The HBChPT contains only terms of the specified order, while both the IR and the EOMS schemes 
contain a series of higher order terms, sometimes termed as recoil corrections. This is illustrated in Fig. |2]for the case of pion- 
nucleon scattering, where we have used filled and half-filled circles to denote the fact that at each chiral order the analytical 
terms of different renormalization scheme for the same loop amplitude are not necessarily the same. 

Furthermore, both the HB and the IR formulations spoil the desired analyticity of loop amplitudes. In addition, the separation 
of a co variant loop result into a singular part and a regular part may make both parts divergent though the sum may be convergent. 
This is the case of the N 3 LO result for the nucleon mass. In Table [I] we have used (y 7 ) to denote this fact. On the other hand, 
although the EOMS seems to be the most desired formulation, it does have its own complications. That is to say, in the EOMS 
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FIG. 2. (Color on-line) Chiral order vs. number of loops for pion-nucleon scattering in BChPT with MS, HB, IR, and EOMS renormalization 
schemes. Red filled circles denote PCB terms. Half-filled circles indicate the fact that the analytical terms of the corresponding BChPT results 
are not necessarily the same as those denoted by filled circles. 



TABLE I. A comparison of different formulations of BChPT. 



power-counting covariance analyticity ultraviolet regularisation 

Ms - V ~ V 

HB sj V 

iR V V (V) 

eoms V V V V 



scheme it is not as straightforward as in the HB and IR approaches to obtain the PCB terms (see, e.g., Refs. ll62ll63ll ). Particularly, 
one has to to make sure that the small variables used in the expansion are defined the same when interpreting or using the LECs 
obtained by others. 

From a practical point of view, the EOMS formulation has recently been shown to converge relatively faster than the IR and 
the HB formulations. Some of these studies will be briefly reviewed in the following section. 



III. RECENT APPLICATIONS IN THE THREE-FLAVOR SECTOR 



In this section, we briefly review some recent applications of the covariant (EOMS) BChPT in the one-baryon sector. In 
particular, we focus on the three-flavor sector of u, d, and s quarks. It should be noted that the EOMS formulation has been 
applied to the baryon-baryon system llr54ll . 



A. Magnetic moments of the octet and decuplet baryons 

1. Virtual octet contributions to the magnetic moments of octet baryons 

In the SU(3) flavor symmetric limit one can relate the magnetic moments of the baryon-octet and the AS transition to those of 
the proton and the neutron. These are the celebrated Coleman-Glashow formulas 16511 . To properly implement SU(3) breaking, 
ChPT should be an appropriate framework to tackle this problem in a systematic fashion. This was first attempted by Caldi and 
Pages even before ChPT as we know today was formulated ll66ll and in the HBChPT framework by Jenkins et al. I67il . It was 
found that the LO SU(3) breaking effects induced by loops are too large and worsen the SU( 3) sy mmetric descriptions. This 
problem has often been used to question the validity of SU(3) BChPT in the one-baryon sector II68h70|1 . 

In the last decade several calculations in HBChPT up to next-to-next-to-leading order (NNLO) have been performed both 
with 1(371 l68l Pnll and without l69tl the inclusion of the baryon decuplet. The large number of LECs appearing at this order 
reduces the predictive power of the theory. The baryon-octet magnetic moments have been calculated using the IR method ll72tl 
and, at NLO, the SU(3)-breaking corrections are still large. Moreover, the agreement with the data is even worse than in 
HBChPT. The size of NLO terms raises the question about the convergence of the chiral series llr3^l - l70ll . 
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TABLE II. Magnetic moments of the octet baryons calculated in LO BChPT, NLO BChPT formulated in the HB, IR and EOMS schemes. All 
the values for the magnetic moments are expressed in units of nuclear magnetons, while b{? and bg are dimensionless. Taken from Ref. (15!:]. 
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2.56 


-1.60 


-0.80 


-0.97 


2.56 


0.80 


-1.60 


-0.97 


1.38 


2.40 


0.77 


0.46 


0(p 3 ) 


HB 


3.01 


-2.62 


-0.42 


-1.35 


2.18 


0.42 


-0.70 


-0.52 


1.68 


4.71 


2.48 


1.01 


IR 


2.08 


-2.74 


-0.64 


-1.13 


2.41 


0.64 


-1.17 


-1.45 


1.89 


4.81 


0.012 


1.86 


EOMS 


2.58 


-2.10 


-0.66 


-1.10 


2.43 


0.66 


-0.95 


-1.27 


1.58 


3.82 


1.20 


0.18 


Exp. 


2.793(0) 


-1.913(0) 


-0.613(4) 


-1.160(25) 


2.458(10) 




-0.651(3) 


-1.250(14) 


± 1.61(8) 










0.2 0.4 0.6 0.8 1 

x 

FIG. 3. (Color on-line) SU(3)-breaking evolution of the minimal \ 2 in the 0(p 3 ) BChPT approaches. The shaded areas are produced by 
varying Mb from 0.8 GeV to 1.1 GeV. This effect lies within the line thickness in the EOMS case, while the HB is insensitive to it. Taken 
from Ref. fBll . 



A study of the magnetic moments of octet baryons in the EOMS formulation was performed in Refs. [|l5l l73ll . Up to NLO, 
there are two LECs appearing in the BChPT calculation, i.e., and b$ . One can fix them by performing a fit of the BChPT 
results to the corresponding experimental data. The results are shown in Table Hfl where x 2 is defined as x 2 — ^2(nth — f^exp) 2 - 
One should note that this is not a proper definition of x 2 , but is used nevertheless to compare with previous studies. Clearly, 
unlike HB and IR, the EOMS results improve the tree level results. On the other hand, the IR results are even more off the data 
than the HB results. 

Furthermore, in Ref. |fl5ll . it was pointed out that the EOMS results show a better convergence behaviour. To see this, one can 
separate the 0(p 2 ) contribution from that of 0(p 3 ): 



li p = 3.47 (1 - 0.257), 
fj, A = -1.27 (1 - 0.482), 
= 3.47 (1 - 0.300), 
fj, s - = -0.93 (1 + 0.025) , 
/i As o = 2.21 (1 - 0.284) . 



H n = -2.55 (1 - 0.175), 
/i s - = -0.93 (1 + 0.187), 
/i s o = 1.27 (1 - 0.482), 
H E o = -2.55 (1 - 0.501) , 



where the first number in the parenthesis indicates the LO contribution and the second number denotes the percentage of the 
NLO contribution relative to that of the LO. It is clear that in the EOMS framework the NLO contribution is at most about 50% 
of the LO, while the same number can be up to 70% in HB or 300% in IR. 

In the chiral limit, all three approaches should give the same description as the tree level results, while in the real world only 
the EOMS results improve the tree level description. Such an evolution is shown in Fig. 2. 



10 



2. Virtual decuplet contributions to the magnetic moments of the octet baryons 

A basic assumption of EFTs is that high-energy dynamics can be integrated out or, in other words, their effects can be 
approximated by contact interactions. In the present case, this means that in the loops virtual baryons except the octet baryons 
can be neglected. In the three-flavor sector, the above assumption may not be valid. The reason is that the lowest-lying decuplet 
is just, on average, 23 1 Me V above the ground-state octet. This energy difference is slightly larger than the pion mass and much 
smaller than the kaon mass. Therefore one needs to check the contributions of the virtual decuplet baryons. Such a work was 
performed in Ref . lf73ll . 

Inclusion of the decuplet baryons is not as straightforward as one would like because the description of higher-spin (s >3/2) 
particles in a relativistic quantum field theory is known to be problematic because of the presence of unphysical lower-spin com- 
ponents. For instance, in the Rarita-Schwinger (RS) formulation lf74ll adopted in this work, the field representation of a massive 
3/2 -particle is a vector-spinor ipu with two unphysical spin- 1/2 components in addition to the spin-3/2 component. In the pres- 
ence of interactions the unphysical degrees of freedom are known to lead to pathologies like non-positive definite commutators 
or acausal propagation for the coupling of the photon IT75l - l77ll . Equivalent problems in phenomenological hadronic interactions 
have also been extensively discussed [78-82]. In the context of ChPT one can use field redefinitions on the conventional chiral 
Lagrangians in order to cast the interactions in a form that is invariant under the transformation ip^ — > ip^ + d^e ll83T[85ll . The 
resulting gauge symmetry ensures to keep active only the physical degrees of freedom HH. Furthermore, there is abundant work 
concerning the inclusion of spin-3/2 resonances in the framework of baryonic effective field theories |[86l - l89ll . 

The baryon-decuplet consists of a SU(3)-flavor multiplet of spin-3/2 resonances that can be represented with the Rarita- 
Schwinger field T M ee T£ de with the following associations: T 111 = A++, T 112 = A+/V3, T 122 = A°/\/3, T 222 = A~, 
T ii3 = Y,*+/VS, T 123 = S*°/V6, T 223 = S*-/\/3, T 133 = S*°/V3, T 233 = H* _ /\/3, and T 333 = Q~. The covariantized 
free Lagrangian is 

C D = ff c {i^ a D a - M Dl ^)T^ bc , (29) 

with M D the chiral limit decuplet-baryon mass and D v Tf c = d v Tf c + {T„) a d T dbc + (T„) b d T^ dc + {T„) c d Tf d . In the last 
and following Lagrangians we sum over any repeated SU(3)-index denoted by latin characters a, b,c, . . ., and (X)% denotes the 
element of row a and column b of the matrix representation of X. 

The conventional lowest-order chiral Lagrangian for the interaction of the decuplet- and octet-baryons with the pseudoscalar 
mesons expanded up to one meson field is 

4b2> = I" £abCf Z de W + Z ^ B c d ^ + hX -' (30) 

where C is the <\>BD coupling, the meson-decay constant and z is an off-shell parameter. An analysis of the constraint 
structure of the interacting theory of Eqs. (|29] [30} yields z = — 1 |p78tl - Nevertheless, the resulting interaction leads to well- 
known problems afflicting the relativistic quantum field theory of 3/2-spinors l!79M8lll . 

The alternative approach of demanding the effective Lagrangians to be spin-3/2-gauge invariant leads, after a field redefinition, 
to the "consistent" (j>BD interaction l8ll[82ll 

C +* D = ^ (9af; " e) ^ + h ' C " (31) 

which is on-shell equivalent to Eq. ( f30l >. In addition, one obtains a second-order <j>cj)BB contact term 

4Sbb = umIf 2 (^{[d^d^l (R^B)}) + (BUd^M], (R^B)]) - 6(Bd^)(d^(R^ B))) , (32) 

where R^ u = i^i^ ua d a + Md^ u and (. . .) denotes the trace in flavor space. The latter Lagrangian is interpreted as carrying 
the spin- 1/2 content of the Lagrangian ( f30l >. This term is eliminated by absorbing it into suitable higher-order LECs. 

Using the EOMS scheme to remove the PCB terms in the loops and the so-called small scale expansion (87] to treat the 
octet and decuplet mass difference, one obtains the results shown in Table Hill Immediately, one notices that in the "consistent" 
decuplet-octet coupling scheme, the inclusion of the virtual decuplet baryons only slightly worsen the octet-only description. 
On the other hand, using the conventional coupling scheme of the decuplet-octet couplings, the inclusion of the virtual decuplet 
deteriorates a lot the description of the magnetic moments of the octet baryons. 

It has been checked that the conclusion does not depend sensitively on the values of the regularisation scale, the decuplet-octet 
mass splitting, and the average baryon mass, as shown in Fig. [4] 
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TABLE III. Baryon octet magnetic moments in BChPT up to 0{p i ). The SU(3)-symmetric description is compared with the different 0(p 3 ) 
BChPT calculations, i.e., the HB and the EOMS results both with (O+D) and without (O) the inclusion of dynamical decuplet baryons. In the 
covariant case the numerical results are obtained using the consistent couplings d31t and the conventional couplings J30b with z = — 1. The 
experimental values are from Ref. l90ll . Taken from Ref. |f73ll . 





Tree level 0{p 2 ) 


Heavy Baryon 0(p A ) 


Covariant EOMS 0{p 3 ) 


Expt. 




O 


O+D 


O 


O+D (conv.) O+D (consist.) 


p 


2.56 


3.01 


3.47 


2.60 


3.18 


2.61 


2.793(0) 


n 


-1.60 


-2.62 


-2.84 


-2.16 


-2.51 


-2.23 


-1.913(0) 


A 


-0.80 


-0.42 


-0.17 


-0.64 


-0.29 


-0.60 


-0.613(4) 


£" 


-0.97 


-1.35 


-1.42 


-1.12 


-1.26 


-1.17 


-1.160(25) 


E+ 


2.56 


2.18 


1.77 


2.41 


1.84 


2.37 


2.458(10) 


E° 


0.80 


0.42 


0.17 


0.64 


0.29 


0.60 






-1.60 


-0.70 


-0.41 


-0.93 


-0.78 


-0.92 


-0.651(3) 


=•0 


-0.97 


-0.52 


-0.56 


-1.23 


-1.05 


-1.22 


-1.250(14) 


AE° 


1.38 


1.68 


1.86 


1.58 


1.88 


1.65 


±1.61(8) 




2.40 


4.71 
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3.92 


5.76 


4.30 




hi 
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2.49 


1.28 
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0.22 
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FIG. 4. Uncertainties of the numerical results of Table ITTTl due to the values of the regularisation scale fi (left panel), the average baryon mass 
Mb (center panel) and the decuplet-octet mass splitting 8 (right panel). The lines represent the \ 2 for the results obtained in HB (dotted line) 
and in the EOMS approach using conventional (dashed line) and consistent (thick solid line) couplings. The grey lines represent the x 2 for the 
case without explicit decuplet resonances in HB (dotted) and in the EOMS formulation (solid). For reference the SU(3)-symmetric description 
(thin solid line) is also shown. Taken from Ref. 17311 



3. Magnetic moments of the decuplet baryons 



The EOMS BChPT has also been applied to study the electromagnetic structure of the decuplet baryons 119 111 . The structure 
of the spin-3/2 particles, as probed by photons, is encoded into four electromagnetic form factors Ir92l : 



(Ttf)\J»\T(p)) = -u a (p'){ [F*(r)^ + !^F*(T)]g a e + [f 3 *(t) 7 " + ^±.FZ{r)] ^t}Mp) 



io-^q v 



(33) 
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TABLE IV. Magnetic dipole moments of the decuplet baryons (in nuclear magnetons) in 0(p 3 ) EOMS BChPT, in comparison with the SU(3)- 
symmetric description and with those obtained in other theoretical appr oache s including the NQM Ir93ll . the RQM 1941. the xQM (HI], the 
XQSM Si, the QCD-SR (H, (extrapolated) 1QCD HHH], large N c fTooh and the HBChPT calculation of Ref. jlOlfl . The experimental 
values are also included for reference ||102|| . Taken from Ref. l9lll . 





A++ 


A 


\0 

i\ 


A — 

A 




Zj 




1 — 1 


1 — 1 


o — 
1 1 


SU(3)-symm. 


4.04 


2.02 





-2.02 


2.02 





-2.02 





-2.02 


-2.02 


NQM [93J 


5.56 


2.73 


-0.09 


-2.92 


3.09 


0.27 


-2.56 


0.63 


-2.2 


-1.84 




4. 10 




u 




i.SZ 


-U.2 / 


-2.J6 


-U.6U 


O A 1 

-z.41 


-2. 3D 


, ./~wa ro^i 
XQM [95] 


C Cf2 


i An 


U 


-J.4 / 


4.1z 


U.J J 


1 C\C 


1 in 
i.iU 


-z.ol 


-Z.lj 


XQSM[96J 


4.85 


2.35 


-0.14 


-2.63 


2.47 


-0.02 


-2.52 


0.09 


-2.40 


-2.29 


QCD-SR [97] 


4.1(1.3) 


2.07(65) 





-2.07(65) 


2.13(82) 


-0.32(15) 


-1.66(73) 


-0.69(29) 


-1.51(52) 


-1.49(45) 


1QCD[98J 


6.09(88) 


3.05(44) 





-3.05(44) 


3.16(40) 


0.329(67) 


-2.50(29) 


0.58(10) 


-2.08(24) 


-1.73(22) 


1QCD[99J 


5.24(18) 


0.97(8) 


-0.035(2) 


-2.98(19) 


1.27(6) 


0.33(5) 


-1.88(4) 


0.16(4) 


-0.62(1) 




large iV c flOOl 


5.9(4) 


2.9(2) 




-2.9(2) 


3.3(2) 


0.3(1) 


-2.8(3) 


0.65(20) 


-2.30(15) 


-1.94 


HBChPT [101] 


4.0(4) 


2.1(2) 


-0.17(4) 


-2.25(19) 


2.0(2) 


-0.07(2) 


-2.2(2) 


0.10(4) 


-2.0(2) 


-1.94 


EOMS BChPT 


6.04(13) 


2.84(2) 


-0.36(9) 


-3.56(20) 


3.07(12) 





-3.07(12) 


0.36(9) 


-2.56(6) 


-2.02 


Expt. [102] 


5.6±1.9 


2.7li; ( 3±1.5±3 
















-2.02±0.05 



Gev(t) 


= (*i(t)- 


^ 2 *(r))H 


- ^tG B2 {t), 




G E 2(t) 


= (*i(t)- 


tF£(t))- 






G M i(t) 


= (^T(r)4 


-f;{t)) + 


4 

7TGM3W, 

5 




G M3 (t) 


= KWi 




I(l + T )(F*(r)- 





where u Q are the Rarita-Schwinger spinors and r = -g 2 /( 4M !>)- One can define the electric monopole and quadrupole, and 
the magnetic dipole and octupole form factors in terms of the F*'s: 

(34) 
(35) 
(36) 
(37) 

At q 2 = 0, the multipole form factors define the static electromagnetic moments, namely, the charge Q, the magnetic dipole 
moment fj,, the electric quadrupole moment Q, and the magnetic octupole moment O 

Q = G E0 (Q) = F?(Q), (38) 

^4 Gm(OH i (Q+ ^ (0)) ' (39) 

Q = il GE2{Q) = i4 (Q - 2^ 3 * (0)) ' (4o) 

° = 2i% G - (0) = 2k ( GMl(0) " ^ + Fm) ) • (41) 

The electromagnetic multipole moments of the spin-3/2 resonances are connected with their spatial electromagnetic distribu- 
tions and, therefore, with their internal structure. Particularly, the EQM and MOM measure the departure from a spherical shape 
of the charge and from a dipole magnetic distribution, respectively. 

Besides the static electromagnetic moments, the slope of the form factors at q 2 = is also of phenomenological interest. In 
particular the one corresponding to Geo is the so-called squared CR: 



4o> = 6- 



dG E0 (q 2 ) _dF*(q 



dq 2 



■2) 



= 6- 



dq 2 



Fm-^2-G E2 {0). (42) 



2M% 2K ' M 2 , 



At NLO, there is only unknown LEC in BChPT, g^, for the magnetic dipole moment, which can be fixed by reproducing the 
magnetic moment of the Q~. Once this is done, BChPT can make a parameter free predictions for all the magnetic moments of 
the other decuplet baryons. The results are shown in Table HVI and compared with those of a number of other studies, including 
the non-relativistic quark model (NQM) JgU, the relativistic quark model (RQM) ll94tl . the chiral quark model ll95ll . the c hiral 
quark soliton m odel ll96tl . the QCD sum rule approach ll97ll . the LQCD approach l98l l99tl . the Large N c approach llOOll . the 
NLO HBChPT llOlll . It is clear that the EOMS BChPT results are in very good agreement with the data. 
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Up to NLO, there is only one unknown LEC g q for the electric quadrupole moments, and one unknown LEC g CI for the charge 
radii of the decuplet baryons. Q If the electric quadrupole moment and the charge radius of a decuplet baryon is known via either 
experiment or LQCD, the unknown LECs can be fixed and one can make predictions for all the other decuplet baryons. The 
results then serve as a stringent test on the validity of covariant BChPT at NLO. This has been explored in Ref. ll9~lll . 



B. Masses and sigma terms of the octet baryons 

1 . Masses of the octet baryons in N 3 LO EOMS BChPT 



Recently, th e lowest-ly ing baryon spectrum, composed of up, down and strange quarks, has been studied by various LQCD 
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collaborations II103M 1 1 111 . Such studies not only improve our understanding of some fundamental questions in physics, such as 



the origin of mass 11 1211 . but also provide a benchmark for future studies based on lattice QCD. In addition, such quark mass, 
volume, lattice spacing dependent simulations also enable us to better understand different aspects of QCD, which in reality has 
fixed quark mass and is continuous and in infinite space time. The latter perspective also allows one to fix the many unknown 
LECs of BChPT, which can not be easily determined, if not possible, by experimental data alone. 

However, because these calculations adopt different lattice setup and all of the m le ad to the same continuum theory, it is 
crucial to test whether the results for baryon masses are consistent with each other 1 1 1 3ll . On the other hand, since lattice QCD 
simulations are performed in a finite hypercube and with larger than physical light-quark masses ll23l . the final results can only 
be obtained by extrapolating to the physical point (chiral extrapolation) and infinite space-time (finite volume corrections). ChPT 
provides a useful framework to perform such extrapolations and to study the induced uncertainties. 

In the past decades, the ground-state (g.s.) octet baryon masses have been studied extensively lfl6ll49ill 14l4l27ll . It is found 
that SU(3) HBChPT converges rather slowly lll28ll . Further more, most calculat i ons are performed only up to NNLO because 
of the many unknown LECs at N 3 LO except those of Refs. llll7lllldll20l[T25l - ll27ll . Regarding ch iral e xtrapolations, Young 
and Thomas JH obtained very good results using the FRR scheme up to NNLO by fitting the LHPC iTTorh and PACS-CS Il05h 
lattice data. In Ref. HI], we applied the NNLO EOMS-BChPT to analyze the same lattice data and found that the EOMS-BChPT 
can provide a better description of lattice data and is more suitable for chiral extra polation purpo ses than HBChPT and NLO 
BChPT. Recently, usin g a p artial summation scheme up to N 3 LO, Semke and Lutz Ill25[ll27l[l29ll found that the BMW Il04ll . 
HSC OH, PACS-CS Eoj, LHPC JTol], and QCDSF-UKQCD lattice results can be well described. 

On the other hand, until very recently, a simultaneous description of all the rif = 2 + 1 lattice data with finite-volume 
effects taken into account self-consistently is still missing. Such a study is necessary for a clarification of the convergence 
problem and for testing the consistency between different lattice simulations. Furthermore, it also provides a good opportunity 
to determine/constrain the many unknown LECs of BChPT at N 3 LO. 

Up to N 3 LO, the Feynman diagrams contributing to the octet baryon masses are shown in Fig. [5] They can be rather easily 
calculated and the detailed formulas can be found in Ref. ||26ll . There are in total 19 LECs that have to be determined, compared 
to only 4 at NNLO. 



/ \ 
/ i 

-• •- 



(c) (d) (e) 

FIG. 5. Feynman diagrams contributing to the octet-baryon masses up to 0(p 4 ) in the EOMS-BChPT. The solid lines correspond to octet- 
baryons and dashed lines refer to NGBs. The black boxes (diamonds) indicate second (fourth) order couplings. The solid dot (circle-cross) 
indicates an insertion from the dimension one (two) meson-baryon Lagrangians. Wave function renormalization diagrams are not explicitly 
shown but included in the calculation. 



The mag netic octupole moments depend on g q and g^. 
5 In Ref. Il27ll . Semke and Lutz showed that their partial summation approach ca n rep roduce the results of the HSC and QCDSF-UKQCD collaborations 
by fitting the BMW, PACS and LHP data. It should be noted, however, that Ref. Il27ll adopted the large jV c argument to fix the values of some LECs and 
neglected the SU(3)-symmetric contributions. 
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FIG. 6. (Color online). Landscape of LQCD simulations of the ground-state octet baryon masses: the PACS-CS (red circles), LHPC (blue 
squares), QCDSF-UKQCD (green diamonds), HSC (yellow upper triangles) and NPLQCD (purple lower triangles) collaborations in the 
2M^ — M% vs. plane (left panel) and in the L vs. M% plane (right panel). The star denotes the physical point with the physical light- and 
strange-quark masses (as implied by leading-order ChPT). Taken from Ref. |26|| . 



The LQCD simulation points are depicted in Fig.|6]in the M 2 vs. 2M\ — M 2 plane and in the M 2 vs. L plane. It is clear 
that to obtain the physical octet baryon masses, one has to extrapolate the LQCD simulation results to the physical light-quark 
masses and infinite space time. On the other hand, as an perturbation theory up to a certain order, the covariant BChPT up to 
N 3 LO is not expected to work up to arbitrarily large light-quark masses and arbitrarily small volumes. Therefore, one has to 
select the LQCD simulations with a limited range of quark masses and volumes. In Ref. |[26i], 1 1 sets of LQCD data are selected 
with Af 2 < 0.5 GeV 2 and M$L > 4 with $ = ir, K, r\. Fitting to the 1 1 sets of LQCD data plus the physical data, one obtains 
the results shown in Table [V] @ An order by order improvement of the description of LQCD mass data is clearly seen from 
the decreasing x 2 /d.o.f for each increasing chiral order. The almost perfect fit at N 3 LO indicates that the LQCD simulation 
results from different collaborations are consistent with each other. Furthermore, the values of the LECs at each order seem to 
be natural, i.e., the values are about 0(1). 

It should be stressed that a self-consistent consideration of finite-volume corrections is necessary to achieve a y 2 /.d. o.f cl ose 
to 1. Without taking into account finite-volume corrections, a x 2 /d.o.f rj 1.9 would have been obtained j26ll . In Ref. lll30ll . it 
is shown that the NNLO EOMS BChPT can describe reasonably well the volume dependence of the NPLQCD data, which is 
consistent with the result of the N 3 LO study performed in Ref. 112611 . 



2. Sigma terms of the octet baryons 

A natural application of the above study is to use the so-obtained LECs to predict the bary on sigma terms, which play important 
roles in understanding chiral symmetry breaking and in direct dark matter searches ill 3 ill - They are related to the light-quark 
mass dependence of the baryon masses via the Feynman-Hellmann theorem, which states: 

dM B 

o-kB = mi(B(p)\uu + dd\B(p)) = mi — , (43) 

ami 

a sB = m s (B(p)\ss\B(p)) = m s ^-, (44) 

am s 

where to; = (m u + md)/2. Using LO meson ChPT, the quark masses can be expressed by the pseudoscalar masses, with 
mi = A/ 2 / (2Bq) and m s = (2Mj^ — A/ 2 )/ (2B n ). Other related quantities, which often appear in the literature, including the 



A second larger set of data, Set-II, has been used to demonstrate the applicability of covariant BChPT in Ref. f26ll . 
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TABLE V. Values of the LECs from the best fit to the LQCD data and the experimental data at 0{p 2 ), 0(p 3 ), and 0(p 4 ) EOMS BChPT. 
Taken from Ref. (53]. 







Set-I 




Set-II 


Fit-C(p 2 ) Fit-C(p 3 ) 


Fitl-C(p 4 ) 


Fit II-C(p 4 ) 


mo [MeV] 


900(6) 


767(6) 


880(22) 


868(12) 


b [GeV _1 l 


-0.273(6) 


-0.886(5) 


-0.609(19) 


-0.714(21) 


b D [GeV -1 ] 


0.0506(17) 0.0482(17) 


0.225(34) 


0.222(20) 


b F [GeV -1 ] 


—0.179(1) 


—0.514(1) 


-0.404(27) 


-0.428(12) 


bi [GeV -1 ] 






0.550(44) 


0.515(132) 


b 2 [GeV -1 ] 






-0.706(99) 


0.148(48) 


63 [GeV -1 ] 






-0.674(115) 


-0.663(155) 


bi [GeV -1 ] 






-0.843(81) 


-0.868(105) 


&5 [GeV -2 ] 






-0.555(144) 


-0.643(246) 


be [GeV -2 ] 






0.160(95) 


-0.268(334) 


67 [GeV -2 ] 






1.98(18) 


0.176(72) 


bx TGeV -2 ! 






473(65") 


—0 0694(4638') 


di [GeV~ 3 ] 






0.0340(143) 


0.0345(134) 


d 2 [GeV -3 ] 






0.296(53) 


0.374(21) 


d 3 [GeV -3 ] 






0.0431(304) 


0.00499(1817) 


di [GeV~ 3 ] 






0.234(67) 


0.267(34) 


d 5 [GeV -3 ] 






-0.328(60) 


-0.445(26) 


d 7 [GeV -3 ] 






-0.0358(269) 


-0.183(12) 


d s [GeV" 3 ] 






-0.107(32) 


-0.307(21) 


X 2 /d.o.f. 


11.8 


8.6 


1.0 


1.6 



strangeness content (ys) and the so-called "dimensionless sigma terms" (Jib, Isb), can be calculated as well: 

= 2{B{p)\ss\B{p)) = mi 2a sB 
(B(p)\uu + dd\B(p)) m s o^b ' 
mi(B(p)\uu + dd\B(p)) u^b 
flB = Mb = Mb ' (46) 

m s {B{p)\ss\B{p)) a sB 

■* sB = U = TT- (47) 

M B M b 

Using the 0(p i ) Fit-I LECs tabulated in Table IVl the pion- and strangeness sigma terms ovs, <t s b for all the octet baryons, 
and the corresponding strangeness content ys, "dimensionless sigma terms" Jib, f S B can be calculated and are shown in Table 
IVII The nucleon pion-sigma term at the physi cal po int, a^N = 43(1) (6) MeV, is in reasonable agreement with the determination 
in the study of the old n — N scattering data Il32ll of a^N = 45 ± 8 MeV, but smaller than the central value of the more recent 
study, (iirN = 59 ± 7 MeV 1 13311 . The ct^n is also i n agr eement with the recent lattice result of the QCDSF collaboration 
( a7rN = 38(12) MeV CH, a^ N = 37(8) (6) MeV HH) and BMW collaboration (a^ N = 39(4) t^ 8 MeV) JH within 
uncertainties, and is cons isten t with the HBChPT result of Ref. rf52"tl . On the other hand, it is larger than the QCDSF resul t 
(a^N = 31(3)(4) MeV) lfl36ll and the N 3 LO BChPT results usin g the partial summation scheme (a^ N = 32(1) MeV) UH, 
but slightly smaller than the JLQCD result (a^N = 50(4.5) MeV) |137ll . 

It should be mentioned that the central value of the nucleon sigma term o^n is smaller than that of Ref. lfl6ll . where the result 
is obtained with the NNLO EOMS-BChPT by fitting the PACS-CS data. Inclusion of virtual decuplet contributions may have 
some non-negligible effects on the predicted baryon sigma terms. At NNLO, it is found that the incl usion of virtual decuplet 
baryons can increase the pion-nucleon sigma term while decrease the strangeness-nucleon sigma term lll38ll . It is interesting to 
check whether such effects still exist at N 3 LO. 

The predicted strangeness-nucleon sigma term is larger than those predicted in Refs. J52lll29ll . It will be interesting to find 
out the origin of such discrepancies by studying the dependence on different formulations of BChPT and on LQCD data. 



7 For a related discussion, see Ref. Il39ll . 
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TABLE VI. The sigma-terms, the strangeness content and the "dimensionless sigma terms" of the octet baryons at the physical point. The first 
error is statistical and the second one is systematic, estimated by taking half the difference between the N' ? LO result and the NNLO result. 
Taken from Ref. jH). 





g^b [MeV] 


g sB [MeV] 


VB 


flB 


fsB 


N 


43(1) (6) 


126(24) (54) 


0.244(47)(110) 


0.0457(11)(64) 


0.134(26)(57) 


A 


19(1)(7) 


269(23) (66) 


1.179(118)(522) 


0.0170(9) (63) 


0.241(21)(59) 


S 


18(2) (6) 


296(21)(50) 


1.369(180)(512) 


0.0151(17)(50) 


0.248(18) (42) 




4(2)(3) 


397(22)(56) 


8.263(4157) (6306) 


0.00303(152)(228) 


0.301(17)(42) 



C. Hyperon vector couplings 

Hyperon semileptonic decays, parameterized by three vector transition form factors (fa, fa, and fa) and three axial form 
factors (#1, g2, an d q^), hav e received renewed interest in recent years due to various reasons. In particular, they prov i de an 
alternative source lfl40t - fl43ll to allow one to extract the Cab ibbo-Kobay ashi-Maskawa (CKM) matrix element V u s Il44lll45ll . 



in addition to kaon semileptonic decays (see, e.g., Ref. 114611 for a recent review), hadronic decays of the r lepton 114711 and the 
ratio T(K + — > ^l + v^)/T{t: + — > /J + ^ M ) lll48ll . The hyperon vector coupling fa (0) plays an essential role in order to extract V us 
accurately. 

Because of the Conservation of Vector Current (CV C) /i(0) is known up to SU(3)-breaking effects, which are of subleading- 
order according to the Ademollo-Gatto theorem lll49ll . The oretical est imates of SU(3 )-bre aking corrections t o fa (0) ha ve been 
performed in various frameworks, including quark models lll50l - ll52ll . large-jV c fits OH, and BChPTM fl53l4l57ll . These 



SU(3)-breaking corrections hav e als o been studied recently in quenched LQCD simulations jl58l Il59 1, and more lately in 
n f = 2 + 1 LQCD simulations lfl60ll . 

In principle, ChPT provides a model independent way to estimate the SU(3)-breaking corrections to /i(0). However, it is 
known that ChPT calculations converge slowly in the three-flavor sector of u, d, and s quarks. This problem becomes even more 
pronounced in the one-baryon sector, where the physics at a certain order can be blurred by the PC restoration procedures, as 
can be clearly seen in the case of the baryon octet magnetic moments lfl5tl . Fortunately, in the case of fa (0), the Ademollo-Gatto 
theorem dictates that up to 0{p A ) no unknown LEC's contribute and, therefore, no PCB terms appear. Consequently, up to this 
order there is no need to apply any PC restoration proce dure s and a BChPT calculation is fully predictive. 

In a recent 0(p 4 ) calculation performed in HBChPT Il56ll . it was shown that the chiral series with only the octet contri butions 
converge slowly while the convergence is completely spoiled by the inclusion of the decuplet ones. In a later work lll57ll . the IR 
BChPT was employed and calculations were performed up to 0(p A ) with only the octet contributions. The slow convergence of 
the chiral series was confirmed but the importance of relativistic corrections was stressed. The first 0(p 4 ) EOMS BChPT study 
was performed in Ref. ll39ll . including the contributions of both virtual octet and virtual decuplet baryons. 

The baryon vector form factors as probed by the charged AS=1 weak current = V us u^^s are defined by 



(B'\V\B) = Vu,S(p') 



^ 2) + jO/^ ( ' 2) + m^«' 2 > 



u(p), (48) 



where q = p' — p. In the SU(3)-symmetric limit, /i(0) is fixed by the conservation of the SU(3)y-charge gy. Furthermore, the 
Ademollo-Gatto theorem states that SU(3)-breaking corrections start at second order in the expansion parameter m s — m 

fa(0)=g v + O((m s -m) 2 ), (49) 



where m s is the strange quark mass and m is the mass of the light quarks. The values of gy are — y |, — —1, y |, 1 for 

A — > p, S° — > p, £~ — > n, S~ — > A, S~ — > S°, and S° — > E + , respectively. In the isospin-symmetric limit only four of these 
channels, which are taken as A — > N, E — > N, 3 — > A, and S — > S, provide independent information. The SU(3)-breaking 
corrections are parameterised order-by-order in the relativistic chiral expansion as follows: 

fa(0) = gy (l + 5^ + 5^ + ...), (50) 

where 5^ and <5^ 3 ) are the LO and NLO SU(3)-breaking corrections induced by loops, corresponding to 0(p 3 ) and 0(p 4 ) chiral 
calculations. 

Up to (D(p A ), one has to calculate the diagrams shown in Figs. [7] and [8] Using the LECs given in Table lVlTl one can calculate 
the contributions of virtual octet baryons and virtual decuplet baryons. The virtual octet contributions are shown in Table [VTlII 
in comparison with the HB and IR BChPT results, when available. It is seen that the HB results are qualitatively different from 
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FIG. 7. Feynman diagrams contributing to the SU(3)-breaking corrections to the hyperon vector coupling /i(0) up to 0(p 4 ). The solid lines 
correspond to baryons and dashed lines to mesons; crosses indicate the coupling of the external current; black dots denote mass splitting 
insertions. The diagrams corresponding to wave function re normalization are not explicitly shown but are taken into account in the calculation. 




, — - , — - , — - /'"*~\ 

-^=^ -J=*#- -H^^U- -1=C^- ^) 

FIG. 8. Feynman diagrams contributing to the LO and NLO SU(3)-breaking corrections to the hyperon vector coupling /i(0), through 
dynamical decuplet baryons. The notations are the same as those of Fig.[7]except that double lines indicate decuplet baryons. 



those of the IR and EOMS results at 0(p i ). In three cases, the sign of 6(2) + 6(3) are different. On the other hand, the EOMS 
and IR results agree much better with each other. Furthermore, one notices that even in the EOMS case, convergence seems 
to be slow. The contributions of virtual decuplet baryons are given in Table HXl in comparison with the HB results. One notes 
that in general the contributions of the virtual decuplet baryons are not negligible and their contributions to the SU(3) breaking 
corrections are positive. 

In Table IXl the complete 0(p 4 ) covariant BChPT predictions are compared with the results of a number of other approaches. 
They seem to agree with those of the chiral quark model and large N c results, but are only partially consistent with the quenched 
LQCD simulations. 

It was pointed out that in Ref. lll60ll that the tendency of the SU(3) breaking correction observed in the n f = 2 + 1 LQCD 
simulations disagrees with predictions of both the latest BChPT result and large N c analysis. It will be of great importance 
to find out how to reconcile the BChPT results with the dynamical LQCD simulations. A careful study of the quark-mass 
dependence of, finite-volume corrections to, and even the finite l attice -spacing dependence of the LQCD simulations may have 
to be performed in order to understand the observed discrepancy lll60ll . 



D. Recent applications in the u, d two-flavor sector 



Although in this short review we mainly focus on the u, d, s three-flavor sector, we must point out that there are many 
interesting and important developments in the two-flavor sector. For instance, the pion-nucleon scattering has been studied up 



TABLE VII. Values for the masses and couplings appearing in the calculation of the SU(3)-breaking corrections to /1 (0). 



D 


0.8 


Mb 


1.151 GeV 


F 


0.46 


M d 


1.382 GeV 


U 


0.0924 GeV 


Mo 


1.197 GeV 


F 


1.17/. 


b D 


-0.0661 GeV" 1 


rrhr 


0.138 GeV 


b F 


0.2087 GeV" 1 


m K 


0.496 GeV 


M D0 


1.216 GeV 


m v 


0.548 GeV 


JM 


0.3236 GeV" 1 


C 


1.0 
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TABLE VIII. Octet contributions to the SU(3)-breaking corrections to /i(0) (in percentage). The central values of the 0(p 4 ) results are 
calculated with /j, = 1 GeV and the uncertainties are obtained by varying fi from 0.7 to 1.3 GeV Taken from Ref. i39ll . 





present work 


HBChPT [156] 


IRChPT [157] 












5 (2) + 5 0) 


5 W +5 (3) 


A - 


+ N 


-3.8 


2+ 12 


-3 6+ 1 ' 2 




2.7 


-5.7 ±2.1 


E - 


+ N 


-0.8 


4 7 +3.8 


o q+3.8 
o.a_ 2 g 




4.1 


2.8 ±0.2 




■» A 


-2.9 


1 7+ 2 - 4 

x - ' -1.8 


-1 2+ 2 - 4 




4.3 


-1.1 ± 1.7 




■» E 


-3.7 


i q+0.3 
J " J — 0.2 


_r r,+0.3 
O.U_ 2 




0.9 


-5.6 ±1.6 



TABLE IX. Decuplet contributions to the SU(3)-breaking corrections to /i(0) (in percentage). The central values of the 0(p 4 ) result are 
calculated with = 1 GeV and the uncertainties are obtained by varying ji from 0.7 GeV to 1.3 GeV. Taken from Ref. |39|| . 







Present work 








HBChPT 






5 (2) 


5 (3) 


S W + 5 (3) 




s m 


J(S) 


5 (2) + 5 (3) 


A - 


> N 


0.7 


3 +01 
°- u -o.i 


q 7+0.1 
J - ' -0.1 




1.8 


1.3 


3.1 


E - 


+ N 


-1.4 




4-8±g : | 




-3.6 


8.8 


5.2 




+ A 


-0.02 


c O+0.4 
J - z -0.3 


c 9+0.4 




-0.05 


4.2 


4.1 




■> E 


0.7 


6.o±i;2 


6 7+ 1 - 9 




1.9 


-0.2 


1.7 



to C(p 3 ) in Refs. Il33l ll6lTl and up to 0(p 4 ) in Ref. |62]. It was shown that the EOMS approach overcomes the limitations 
that previous chiral analyses of the 7rA^ scattering amplitude had, providing an accurate description of the partial wave phase 
shifts of the Karlsruhe-Helsinki and George-Washington groups up to energies just below the resonance region. Further more , 
the EOMS BChPT is shown to exhibit a good convergence J62fl . The prot on Co mpton scattering has been studied in Ref. lll62ll 
and the nucleon and A(123 2)-is obar electromagnetic form factors in Ref. lll63ll . The nucleon to A axial transition form factors 
have been studied in Ref. Il64ll . The nucleon mass has been studied up to 0(p 6 ) l55l l56ll . the axial, induced pseudoscalar, 
and pion-nucleon form factors up to 0(p 4 ) 15711 . and the electromagnetic form factors of the nucleon up to C(p 4 ) 1580 . The 
amplitude for ordinary muon capture on the proton has been studied as well l63ll . 

S ome of these studies showed clear improvement of the EOMS scheme over the IR or HB formulations (see, e.g., Refs. d 1 3 31. 
Il63l0 . while others did not (see, e.g., Ref. ll63lo . Nevertheless, it is always preferable to perform a study with a fully covariant 
ChPT supplemented with a PC restoration procedure that does not spoil the analyticity of the loop results, e.g., the EOMS 
BChPT. 



E. Extensions to the heavy-light system 

Given the formal similarity between the heavy-light mesons and the ground-state baryons, in particular, regarding chiral 
descriptions of their static properties and scatteri ng w ith the NGBs, the covariant formulation of ChPT has bee n ex tended to 
study the decay constants of D and B mesons lfl7lll65ll . and the scattering lengths of the NGBs off the D mesons Il66ll . Indeed, 
the covariant formulation also shows improved convergence in comparison with the conventional heavy-meson (HM) ChPT. 
From Fig. [9] it is clear that with the same number of LECs, the covariant ChPT describes better the light-quark mass dependence 



TABLE X. SU(3)-breaking corrections (in percentage) to /i(0) obtained in different approaches. Taken from Ref. l39ll . 





Present work 


Large 


Quark model 


quenched LQCD 






Ref. [142] 


Ref. [150] 


Ref. [151] 


Ref. [152] 




A - 


> N 


i +1 - 3 

U - A -1.0 


2±2 


-1.3 


-2.4 


0.1 




E - 


+ N 




4±3 


-1.3 


-2.4 


0.9 


-1.2 ± 2.9 ± 4.0 [158] 




■+ A 


4.0t 2 ;? 


4±4 


-1.3 


-2.4 


2.2 






+ E 


1 7+2.2 
l - '-1.6 


8±5 


-1.3 


-2.4 


4.2 


-1.3 ±1.9 [159] 
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0.0 0.2 0.4 0.6 0.8 1.0 

m q /m s 

FIG. 9. (Color online) Light-quark mass dependence of f D (solid lines) and f Ds (dashed lines). The black lines show the results of the 
NLO HMxPT and the blue lines the results of the covariant NLO ChPT. The red lines are the continuum extrapolations of the HPQCD 
collaboration lll68ll . The ratio r = m q /m s is related to the pseudoscalar meson masses at leading chiral order through m% = 2Bom a r and 
m 2 K — Bom B (r + 1) with Bo = rn^/(m u + ma) and m q = (m u + md)/2, where m u , mj, and m s are the physical up, down, and strange 
quark mass. Taken from Ref. lfl7tl . 



of the LQCD results than the non-relativistic HM ChPT. 

Recently, the IR f ormulation of ChPT has been applied to study the pseudoscalar meson and heavy vector meson scattering 
lengths as well lll67h . 



IV. SUMMARY AND OUTLOOK 

In this work, we have briefly reviewed some recent developments in covariant baryon chiral perturbation theory (BChPT). 
After a short introduction to the power-counting-breaking problem appearing in the formulation of BChPT and the existing 
recipes to recover a systematic power-counting, we have summarized the advantages and limitations of each method and have 
shown a number of successful applications of the EOMS scheme in the u, d, and s three-flavor sector, including the magnetic 
moments of the octet and decuplet baryons, the masses and sigma terms of the octet baryons, and the hyperon vector couplings. 
The developments in the u and d two-flavor sector, as well as the recent extensions to the heavy-light system are also briefly 
covered. 

It is shown that the covariant BChPT description of the magnetic moments of the octet and decuplet baryons and the masses of 
the octet baryons are very satisfactory, but for the strangeness-nucleon sigma term and the hyperon vector couplings, more works 
are still needed to reconcile the BChPT predictions with the LQCD data or to understand the large discrepancy between different 
studies within the BChPT framework. In this regard, effects of higher-energy degrees of freedom have to be examined more 
carefully, e.g., the decuplet baryons. Although virtual decuplet contributions are found to be small in the case of the magnetic 
moments of the octet baryons, their contributions to the hyperon vector couplings are not negligible. As having been shown 
in the HBChPT, one may have to study the effects of virtual decuplet baryons case by case. Furthermore, the extension to the 
heavy-light system should be further explored. 

Compared to its non-relativistic and infrared counterparts, the EOMS formulation seems to be more suitable for the three- 
flavor sector and for descriptions of LQCD simulation data. One should note, however, that a systematic study of different 
physical observables in the three-flavor sector using the EOMS formulation is still missing. Such a study is necessary not only 
for cross-checking the validity of the EOMS BChPT but also for truly testing the predictive power of the covariant BChPT. In 
view of the rapid progress being made in both experiment and LQCD simulations, such studies are strongly encouraged. 

The manifestly Lorentz covariant formulation of BChPT can also be applied to baryon-baryon scattering. The investigations 
in this direction could further deepen our understanding of the nucleon-nucleon, hyperon-nucleon, and hyperon-hyperon forces. 
With the relevant low-energy constants determined either with the help of LQCD simulations or through fitting experimental 
data, the baryon-baryon forces can be used in studies of the properties of atomic nuclei/hypernuclei. 
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